Commentary on rainbow-ladder truncation for excited states and exotics 
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Ground-state, radially-excited and exotic scalar-, vector- and flavoured-pseudoscalar- mesons are 
studied in rainbow-ladder truncation using an interaction kernel that is consonant with modern 
DSE- and lattice-QCD results. The inability of this truncation to provide realistic predictions for 
the masses of excited- and exotic-states is confirmed and explained. On the other hand, its appli- 
cation does provide information that is potentially useful in proceeding beyond this leading-order 
truncation, e.g.: assisting with development of projection techniques that ease the computation of 
excited state properties; placing qualitative constraints on the long-range behaviour of the inter- 
action kernel; and highlighting and illustrating some features of hadron observables that do not 
depend on details of the dynamics. 
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I. INTRODUCTION 

Meson spectroscopy is a keystone of extant and forth- 
coming programmes at numerous facilities worldwide, 
e.g.: the Beijing Spectrometer; the COMPASS detec- 
tor at CERN; Hall-D at Jefferson Laboratory; the Japan 
proton accelerator research complex (J-PARC); and the 
PANDA detector at GSI. Each identifies an essentially 
identical primary motivation; namely, seeking answers to 
two fundamental questions within the Standard Model: 
What matter is possible; and How is it constituted? 
The subtext is quantum chromodynamics (QCD), the 
strongly-interacting part of the Standard Model, and the 
unique nature of the forces it seems to produce. With 
QCD, Nature has prepared the sole known example of a 
strongly-interacting quantum held theory that is defined 
by degrees-of-freedom which cannot directly be detected; 
i.e., they are confined. One of the greatest challenges in 
modern physics is to comprehend and explain the phe- 
nomenon of confinement. 

Following Ref. Q, confinement in mesons has typically 
been associated with a linearly rising potential between 
the quark-antiquark pair Q . There are sound reasons for 
using such potential model phenomenology in the study 
of heavy quarkonia However, that is not true for 
light-quark systems. The static potential measured in 
simulations of lattice-QCD is not related in any known 
way to the question of light-quark confinement. Light- 
quark creation and annihilation effects are fundamentally 
nonperturbativc. Hence it is impossible in principle to 
compute a potential between two light quarks 0, [1[ . On 
the other hand, confinement can be related to the ana- 
lytic properties of QCD's Schwinger functions so 
the question of light-quark confinement may be trans- 
lated into the challenge of charting the infrared behavior 
of QCD's /3-function. 

To a large degree, this is also true of explaining dynam- 
ical chiral symmetry breaking (DCSB), a phenomenon 



which has an enormous impact on the measurable prop- 
erties of mesons and baryons [T2, EH ■ It is known that 
DCSB; namely, the gen eration of mass from nothing, does 
occur in QCD [15l - [l7l ]. It arises primarily because a 
dense cloud of gluons comes to clothe a low-momentum 
quark [ll|, [l8| . This is readily seen by solving the Dyson- 
Schwinger equation (DSE) for the dressed-quark propa- 
gator; i.e., the gap equation. However, the origin of the 
interaction strength at infrared momenta, which guaran- 
tees DCSB through the gap equation, is currently un- 
known. This relationship ties confinement to DCSB. 
The crucial role of DCSB means that reliable informa- 
tion about the /3-function can only be obtained via a 
symmetry-preserving treatment of the bound-state prob- 
lem that is capable of veraciously expressing DCSB. The 
DSEs provide such a framework [7H13I] and will be em- 
ployed herein. 

A considerable body of recent work (e.g., Refs. [ll|, 
[H, Il9r- 30] ) has shown that in order to gain sensitiv- 
ity to the long-range part of the interaction, one should 
minimally study the properties of mesons with signifi- 
cant rest-frame quark orbital angular momentum, such 
as scalar- and pseudovector-mesons, the radial excita- 
tions of pseudoscalar- and vector-mesons, and tensor 
mesons. A challenging aspect of this problem is that the 
leading-order (rainbow-ladder) in the most widely used 
symmetry-preserving DSE truncation scheme pH |32| 
fails to adequately express the full power of DCSB in 
the kernels of the bound-state Bethe-Salpeter equations 
(BSEs) [H, HI H3. Consequently, the results pro- 
duced for systems other than ground-state flavoured- 
pseudoscalar- and vector-mesons have most often been 
qualitatively and quantitatively incorrect. 

Is there any reason then to revisit the problem of the 
spectrum of excited and exotic mesons using the rainbow- 
ladder truncation? The answer is "no," if the goal is 
to extract quantitatively reliable information about the 
infrared behaviour of QCD's /3-function. On the other 
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hand, the answer is "yes," if one can exploit the trunca- 
tion's simplicity in order to: identify features of excited 
and exotic states that are plausibly independent of the 
truncation; or techniques that can be useful in connec- 
tion with more sophisticated truncations. Such is our 
aim herein. 

In Sec.|TT] we present the gap- and Bethe-Salpeter- 
cquations in the symmetry-preserving rainbow-ladder 
truncation, explain the structure of their solutions and 
define their kernels. Scction lnTl rcDorts and interprets our 
numerical results, which include: masses and decay con- 
stants; an investigation of the relative importance of var- 
ious Dirac structures within meson Bethe-Salpeter am- 
plitudes; and an exploration of the pointwise behaviour 
and sign of the leading invariant amplitudes. Section ITVl 
is an epilogue. 



confinement and/or DCSB are realised in solutions of the 
gap equation. 

The solution of the gap equation is a dressed-quark 
propagator 



S(p) 



1 



ij ■ pA(p 2 ,( 2 ) + B(p 2 X 2 ) il -P + M(p 2 ) 



(4) 

which is obtained from Eq. (TTJ) augmented by a rcnormal- 
isation condition. A mass-independent scheme is a useful 
choice and can be implemented by fixing all renormalisa- 
tion constants in the chiral limit. Notably, the mass func- 
tion, M(p 2 ) = B(p 2 ,( 2 )/A(p 2 ,( 2 ), is independent of the 
renormalisation point, (; and the renormalised current- 
quark mass is given by 



m< = Z m (C,A)m bm (A) 



bm 



(5) 



II. GAP AND BETHE-SALPETER EQUATIONS 

The renormalised rainbow-gap- and ladder-Bethe- 
Salpeter-equations are, respectively: 

Sip)- 1 = Z 2 (i 1 -p + m hm ) 

+ z 2 l%(i)e 2 D^(i)^s(p - (i) 

,A 

T M (k; P) = -Zl / Q{{k - q) 2 ) {k - q) 2 D l ™ c (k ~ q) 
Jq 

X^S(q + )T M (q;P)S(q_)^, (2) 

where: we use a Euclidean metric (T2J; J £ A := J A j0^ 
represents a Poincare-invariant regularization of the in- 
tegral, with A the ultraviolet regularization mass-scale; 
Z2 (C, A) is the quark wavefunction renormalisation con- 
stant, whose location and strength in these equations 
may be understood from Refs. [12, HH; D l ™ e (£) is the 
Landau-gauge free-gauge-boson propagator^ one can 
choose q± = q ± P/2 without loss of generality in this 
Poincare covariant approach; and 



£ 2 g(£ 2 ) = £ 2 g m (f) + 47ra pQC D(0 



(3) 



specifies the interaction, with a P QCD(fc 2 ) a bounded, 
monotonically-decrcasing regular continuation of the 
perturbative-QCD running coupling to all values of 
spacelike- £ 2 , and Gm(^ 2 ) an Ansatz for the interaction 
at infrared momenta, such that Gib,(£ 2 ) <§; <S p qcd(^ 2 ) 
W 2 > 2GeV 2 . The form of G m (£ 2 ) determines whether 



1 Landau gauge is used for many reasons I351I36H . for example, it is: 
a fixed point of the renormalisation group; that gauge for which 
sensitivity to model-dependent differences between Ansatze for 
the fermion-gauge-boson vertex are least noticeable; and a co- 
variant gauge, which is readily im plem ented in numerical simu- 
lations of lattice regularised QCD |37l . 



wherein Z4 is the renormalisation constant associated 
with the Lagrangian's mass-term. Like the running cou- 
pling constant, this "running mass" is a familiar con- 
cept. However, it is not commonly appreciated that 
is simply the dressed-quark mass function evaluated at 
one particular deep spacclikc point; viz, 



M(C 2 ) 



(6) 



The renormalisation-group invariant current-quark 
mass may be inferred via 



rhf = lim 

p 2 — >oc 



2 A 2 



P 



QCD 



M f (p 2 ) , 



(7) 



where / specifies the quark's flavour, j m = 12/(33 — 
2Nf a ): Nf a is the number of quark flavours em- 
ployed in computing the running coupling; and Aqcd 
is QCD's dynamically-generated renormalisation-group- 
invariant mass-scale. The chiral limit is expressed by 



0. 



Moreover, 



VC 2 »A 2 3CD , 



M fl ( P 2 =( 2 ) _ "fr 
m h 



M h (p 2 =Q 2 ) 



(8) 



(9) 



We would like to emphasise, however, that in the pres- 
ence of DCSB the ratio M f x {p 2 ) / 'M f 2 (p 2 ) is not indepen- 
dent of p 2 : in the infrared; i.e., Vp 2 < Aq CD , it then 
expresses a ratio of constituent-like quark masses, which, 
for light quarks, are two orders-of-magnitudc larger than 
their current-masses and nonlinearly related to them 
[H, HI]- (See, e.g., the discussion following Eq. ([T5]).) 

The BSE is an eigenvalue problem for the meson 
masses-squared; i.e., in a given channel Eq. ([2J has so- 
lutions only at particular, isolated values of P 2 = —m 2 M . 
At these values, solving the equation produces the as- 
sociated meson's Bcthc-Salpctcr amplitude, which can 
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then be used in the computation of observable proper- 
ties. Herein we consideiQ flavoured-pseudoscalar-, scalar- 
and vector- meson ground-, radially-excited- and exotic- 
states, so that the following amplitudes arise: 

4 

r> =0 - (*; p) = Y1 ^ T o- ( fc > p ) K- p ), (io) 

i=l 
4 

r 0+ (t, p) = J2 T o+ ( fc > P) F o+ n (ii) 

i=l 
8 

r 1 -(fc;P)=^r^(fc,P)^_(fc;P), (12) 

i=l 

with (oj := a M - P p a • P/P 2 ) 



confinement and DCSB. Moreover, it is consistent with 
modern DSE and lattice studies, which indicate that the 
gluon propagator is a bounded, regular function of space- 
like momenta that achieves its maximum value on this 
domain at s = [42H4^ |. and the dressed-quark-gluon 
vertex does not possess any structure which can qualita- 
tively alter this behaviour [4^, |4(| . Notably, as illustrated 
in Ref. [30j . the parameters D and u> are not independent: 
with Dlo = constant, one can expect computed obscrv- 
ables to be practically insensitive to uj on the domain 
u G [0.4, 0.6] GcV. 



'o- 



'o+ 



ii 



D ■ 



7-P, r 2 + =fc.Pr 2 _ 



k- Pt, 



0+> 'o+ 



'0+ 
■ T 



7?_ =*fcJfc-P 7 .P, 
n 4 - =*[7j7-^7-fc T + fcj7-n 



fc • P[ 7 J 7 T ■ k - 7 ■ fc T 7 J] 
(fc T ) 2 (7j7-^-7-^7j) 



-2fcJ 7 • fc T 7 • P, 



e 7 -F 7 -P. 



(13a) 
(13b) 



P 2 7 • fc - fc ■ P 7 ■ P, (13c) 
(13d) 
(13e) 
(13f) 
(13g) 
(13h) 
(13i) 
(13j) 



(13k) 
(131) 



The canonical normalisation condition (see, e.g., Eq. (27) 
in Ref. [l^] or, more generally, Ref. |4l|) constrains the 
bound-state to produce a pole with unit residue in the 
quark-antiquark scattering matrix. 

It remains only to specify the interaction in order to 
proceed. We use that explained in Ref. [3(J; viz., 



5(s) 



8tt 5 



De~ s/u 



37r 2 7m T(a) 



ln[r + (l + S /A 2 QCD ) 2 ] 



(14) 



where: 7m = 12/25, A QCD = 0.234 GeV; r = e 2 - 1; 
and P(s) = {1 - cxp(-s/[4m 2 ])}/s, m t = 0.5 GcV. This 
interaction preserves the one-loop renormalisation-group 
behavior of QCD in the gap- and Bethe-Salpeter equa- 
tions [lH , and the infrared behaviour can serve to ensure 



III. NUMERICAL RESULTS FOR 
BOUND-STATES PROPERTIES 



A. Ground states 



Using the method of Ref. [47J , we solved the gap equa- 
tion for light u = d quarks and the s-quark, with 
their current-quark masses fixed by requiring that the 
pion and kaon BSEs produce rj 0.138 GeV and 
rriK ~ 0.496 GeV. This is straightforward in rainbow- 
ladder truncation because there is no coupling between 
the separate gap equations and no feedback from the 
BSEs Hil; and yields 



3.4 MeV, = 82MeV 



(15) 



quoted at our renormalisation point £ = 19 GeV, a value 
chosen to match the bulk of extant studies. These values 
correspond to rcnormalisation-group-invariant masses of 
inu,d = 6 MeV, rh s = 146 MeV, one-loop-evolved masses 
of ml2f = 5 MeV, mJ GcV = 129 MeV; and give 
m s / m u = 24. They are consequently comparable with 
contemporary estimates by other means [49[. NB. With 
ui = 0.6 GeV, Mf/Mf = 1.52 < m s /m u , where the 
constituent-quark mass Mj := {s\s > 0, s = Mj(s)}. 

In Table U we report selected results related to ground- 
state pseudoscalar-, scalar- and vector- mesons. The me- 
son masses are obtained in solving the BSEs. Regarding 
the other meson quantities, in terms of the canonically 
normalised Bcthc-Salpetcr amplitudes and with 



2 Masses and other properties of charge-neutral pseudoscalar 
mesons arc affected by the non-Abelian anomaly. In the BSE 
context, this is discussed in Ref. l4Cfl . Since the non-Abelian 
anomaly is a correction to rainbow-ladder truncation that is qual- 
itatively different to the focus of our study, herein we specialise 
to flavoured pseudoscalars. 



Xjf2 (fc;P) = S fl (k + )T JP (k;K)S f2 (k.), (16) 
where fi, ji arc the meson's valence-quark and 
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TABLE II. Masses obtained with Eq. (0, Dui = 
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TABLE I. Results obtained using the interaction in Eq. (|14[) 
with Dlj = (0.8 GeV) 3 . The current-quark masses at £ = 
19 GeV are given in Eq. (|15|) . Dimensioned quantities are re- 
ported in GeV. For comparison, some experimental values 
are [H]: f w = 0.092 GeV, m n = 0.138 GeV; f K = 0.113 GeV, 
m K = 0.496 GeV; f p = 0.153 GeV, m p = 0.777 GeV; and 
U = 0.168 GeV, m = 1.02 GeV. NB. The scalar mesons 
listed here are not directly comparable with the lightest 
scalars in the hadron spectrum because the rainbow-ladder 
truncation is a priori known to be a poor app roximation in 
this channel: nonresonant corrections [25l . [2811 and resonant 
final-state interactions are both important [391 ]. 



(22) 
(23) 



-antiquark, respectively, one has [19|, |5Cj, [51| 

/o- P m = Z 2 tr CD J *757mX 0i - 2 ( fc ; P ) 

ipl- = Z A tr CD / 75X -( fc ; P )» 
12 Jk 12 

/o+ P p- = ^2 tr CD J HnX + 2 (k; P) , 
,A 

Pn+ = -^4tr C D/ X +(k\P) , 

^r 2 m ir2 = i Z2 trcD / TfXir a p ) • 



(17) 
(18) 
(19) 
(20) 
(21) 



Furthermore, the products f Q ± p Q ± describe in-meson 
condensates [l^, HO, H3] - 

B. Radial excitations and exotics 

In addition to properties of the ground-states, we have 
computed selected quantities associated with J = 0, 1 
radial excitations and exotics. In the Poincare covariant 
DSE treatment, exotic states appear as poles in vertices 
generated by interpolating fields with "unnatural time- 
parity" ■ Results are presented in Table HU The last 
column in the Table was prepared as follows. We fitted 
the entries in each row to both m(w) = constant and 



m(u>) = uj(cq + ciw), 



(24) 



then computed the standard-deviation of the relative er- 
ror in each fit, eg for the constant and a-i for Eq. (|24[) . 
and finally formed the ratio: 020 = 02 /co- 
in preparing the table we used Doj = (1.1 GeV) 3 . This 
has the effect of inflating the ir- and p-meson ground- 
state masses to a point wherefrom corrections to rainbow- 
ladder truncation can plausibly return them to the ob- 
served values [10, [57}. It is therefore notable that, in 
contrast to Table U the value reported for m a in Table UT1 



The Table confirms that, with Dlj = constant, observable 
properties of ground-state scalar-, vector- and flavoured- 
pscudoscalar-mcsons computed with Eq. (|14[) arc practi- 
cally insensitive to variations of w £ [0.4, 0.6] GeV. 

It is noteworthy, and readily verified using entries 
in the Table, that the pseudoscalar- and scalar- meson 
masses satisfy the following identities, exact in QCD 



3 Notwithstanding complexities associated with the structure of 
light-quark scalars l39l l52t |53|I , the identity written here applies 
to any scalar meson that can be produced via e+e — annihilation. 
It is not of experimental significance, however, if the pole is deep 
in the complex plane. 
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FIG. 1. Pseudoscalar mesons. Relative difference between 
the mass computed with all the amplitudes in Eq. (|10|l and 
that obtained when the identified i > 2 amplitude is omitted: 
circles - ground-state pion; squares - J PC = CP - exotic; 
and diamonds - first pseudoscalar radial excitation. In all 
cases, uj — 0.6 GeV, Duj = (1.1 GeV) 3 . There is only minor 
quantitative variation with uj £ [0.4, 0.6] GeV. NB. The i = 1 
amplitude is never omitted, it specifies the reference value. 
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FIG. 2. Scalar mesons. Relative difference between the mass 
computed with all the amplitudes in Eq. (|ll|l and that ob- 
tained when the identified i > 2 amplitude is omitted: circles 
- ground-state u = d scalar; squares - J FC = + ~ exotic; 
and diamonds - first pseudoscalar radial excitation. In all 
cases, uj = 0.6 GeV, Duj = (1.1 GeV) 3 . There is only minor 
quantitative variation with uj £ [0.4, 0.6] GeV. NB. The i = 1 
amplitude is never omitted, it specifies the reference value. 



matches estimates for the mass of the dressed-quark-core 
component of the cr-meson obtained using unitarised chi- 
ral perturbation theory [H, H3] ■ 

A comparison between the w-dependence of ground- 
state properties and those of excited- and exotic-states 
was drawn in Ref. [30l | and we only summarise it here. 
Ground-state masses of light-quark pseudoscalar- and 
vector-mesons are quite insensitive to uj G [0.4, 0.6] GeV. 
Any minor variation is described by a decreasing func- 
tion. In the case of exotics and radial excitations, the 
variation with uj is described by an increasing function 
and the variation is usually significant. This is readily 
understood. The quantity r u := 1/uj is a length-scale 
that measures the range over which the infrared part of 
Eq. (|3|), Gm , is active. For uj = this range is infinite, but 
it decreases with increasing uj. One expects exotic- and 
excited-states to be more sensitive to long-range features 
of the interaction than ground-states and, additionally, 
that their masses should increase if the magnitude and 
range of the strong piece of the interaction is reduced 
because there is less binding energy. 

Table |TT] confirms a known fault with the rainbow- 
ladder truncation; viz., whilst it binds in exotic channels, 
it produces masses that arc too light, just as it does for 
axial-vector mesons. It is similarly noticeable that m ni 
is far more sensitive to variations in uj than is m pi ; and 
although m vi < m pi for u = 0.4 GeV, the ordering is 
rapidly reversed. Thus, in conflict with experiment, one 
usually finds m 7Tl > m Pl in rainbow-ladder truncation. 
This, too, is a property of the truncation, which is insen- 
sitive to the details of G(k 2 ); e.g., the same ordering is 
obtained with a momentum-independent interaction [571 ] . 



C. Structure of bound states 

In order to develop insight, both into the structure of 
excited- and exotic-states, and for progressing beyond 
rainbow-ladder truncation, it is useful to know which 
of the invariant amplitudes in Eqs. ([TU)) - (|12[) are domi- 
nant. One useful measure of an amplitude's importance 
is the contribution it makes to a given meson's mass. 
Figure [T] displays the result for pseudoscalar mesons: in 
all cases a good approximation is obtained by retaining 
Fj- and . This outcome is in agreement with ex- 
tant ground-state computations but extends those 
rainbow-ladder conclusions to excited- and exotic-states. 
Evidently, there is little here to distinguish between the 
exotic and the radial excitation. Curiously, plays 
a role of similar magnitude in each state and the am- 
plitudes F 3 _ and are always largely unimportant. 
These last two, in this instance small, amplitudes are 
those most directly associated with nonzero quark orbital 
angular momentum in the meson's rest-frame. 

For scalar mesons, on the other hand, one reads from 
Fig. [2] that FJ+, and Fq + should be included if a re- 
liable approximation is to be obtained. The latter two 
amplitudes are directly associated with significant rest- 
frame quark orbital angular momentum. Notably, in 
quantum mechanical models, scalar mesons are identified 
as 3 Po states, in contrast to ^So for pseudoscalar mesons. 

The vector meson ( 3 Sj) situation is displayed in Fig. [3] 
In agreement with Ref. [58| , a good approximation for the 
vector-meson ground-state is obtained by retaining FjL , 
, F x 5 _ . The last two amplitudes are associated with 
P-wave components in the rest-frame. However, for the 
first radial excitation, Ff_ is also important: this am- 
plitude is directly associated with a P-wave component 
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FIG. 3. Vector mesons. Relative difference between the mass 
computed with all the amplitudes in Eq. (|12[l and that ob- 
tained when the identified i > 2 amplitude is omitted: circles 
- ground-state u = d vector; squares - J PC = 1 h exotic; 
and diamonds - first vector radial excitation. In all cases, 
uj = 0.6 GeV, Duo = (1.1 GeV) 3 . Whilst there are quantita- 
tive changes with uj, the pattern of amplitude importance is 
unchanged. NB. The i = 1 amplitude is never omitted, it 
specifies the reference value. 




in the radially-excited vector-meson's rest frame. These 
observations suggest that a BSE might be built which 
projects selectively onto the first radially excited state. 

The additional information contained in these figures 
indicates that the shortcomings identified above, of the 
rainbow-ladder truncation for states other than ground- 
state vector- and flavoured-pseudoscalar-mesons, can be 
attributed to this truncation's inadequate expression in 
the Bcthc-Salpcter kernels of effects which in quantum 
mechanics would be described as spin-orbit interactions. 
Namely, treating the quark-gluon vertex as effectively 
bare in both the gap- and Bethe-Salpeter-equations leads 
to omission of critically important helicity-flipping inter- 
actions that are dramatically enhanced by DCSB, as dis- 
cussed in Rcfs. m, m, nu. 

One may readily expand on this. For example, vec- 
tor meson bound states possess nonzero magnetic- and 
quadrupole-moments [59[. This fact, Fig. [3] and the as- 
sociated discussion together indicate that there is ap- 
preciably more dresscd-quark orbital angular momen- 
tum within these states than within pseudoscalar mesons. 
Hence, spin-orbit repulsion could significantly boost m Pl 
and thereby produce the correct level ordering; viz., 
m pi > m Vl . Moreover, since exotic states appear as 
poles in vertices generated by interpolating fields with 
"unnatural time-parity," the importance of orbital an- 
gular momentum within these states is magnified. These 
comments apply with equal force to tensor mesons, which 
cannot be formed without rest-frame quark orbital angu- 
lar momentum. 




- 



0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

k 2 [GeV 2 ] k 2 [GeV 2 ] 

FIG. 4. Pseudoscalar mesons, ^-dependence of low-order 
Chebyshev-projections of leading invariant amplitude for 
ground-, radially-excited- and exotic-states: upper four pan- 
els, ground and radial; lower four panels, ground and exotic. 
In all panels, solid - zeroth moment, ground-state; dashed 
- leading moment, comparison state; dash-dot - sublead- 
ing moment, comparison state. Row-1, left, uj = 0.4 GeV; 
Row-1, right, uj = 0.5 GeV; Row-2, left, uj — 0.6 GeV; and 
Row 2, right, uj = 0.7 GeV. This pattern is repeated in 
the next two rows. The normalisation is chosen such that 
^(p 2 = 0) = 1; and Duj = (1.1 GeV) 3 . 



proact0 is provided by the essentially nonperturbative 
DSE truncation scheme whose use is illustrated most 
fully in Ref. [28j . That symmetry-preserving scheme 
deeply embeds effects associated with DCSB into the 
Bcthe-Salpeter kernel. 



At present the best hope for a realistic description 
of the meson spectrum within a Poincare covariant ap- 



4 A lattice-QCD perspective on the meson spectrum may be drawn 
from Ref. loot. 
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D. Connecting amplitudes with observables 

Whilst not directly observable, the momentum- 
dependence of meson Bethe-Salpeter amplitudes is a cru- 
cial determinative factor in the computation of measur- 
able quantities. In Figs. 0] and [SJ therefore, we depict the 
w-dependence of a few low-order Chebyshev moments of 
the leading invariant amplitude for the pscudoscalar and 
vector mesons: 



n F M (p 2 ):=- / dxy/l-x 2 U n (x)F M (k 2 ,x;P 2 ), 



(25) 

where k ■ P = xVk 2 P 2 and U n {x) is a Chebyshev poly- 
nomial of the second kind. NB. For pscudoscalar and 
vector states with natural C-parity, only the even mo- 
ments are nonzero, whereas it is the odd moments which 
are nonzero for the exotic partners of these states. 

The upper four panels in Fig. 2] compare the ampli- 
tudes of the ground-state and first-radially-excited pseu- 
doscalar mesons. The ground-state is clearly insensitive 
to W. However, as hoped for and anticipated, the ra- 
dial excitation reacts strongly to variations in lu. Most 
notable is the suppression of °E ni with decreasing lu, to 
be replaced by an increasingly large E V1 . Indeed, at 
u) = 0.4 GeV, °E m is almost negligible and possesses two 
zeros, instead of the single zero expected in the amplitude 
of a first radial excitation since the work of Ref. [20[ • ^ n 
such circumstances, the radial excitation may even pos- 
sess a smaller charge radius than the ground state [21 1. 

In our view these features signal that values of lu < 
0.5 GeV in Eq. (fT4|) are unphysical; i.e., the long-range 
behaviour of a realistic /3-function cannot dramatically 
suppress the radial excitation's leading amplitude nor 
induce it to have a second zero. This perspective is 
supported by the following considerations. Neither the 
homogeneous BSE nor the canonical normalisation con- 
dition fix the sign of the Bethe-Salpeter amplitude at 
k 2 = 0. As in quantum mechanics, this is arbitrary and 
cannot affect observables. Another parallel with quan- 
tum mechanics is also relevant. Namely, for a ground- 
state, the sign of the radial wave function at the origin 
in configuration space is the same as that of its analogue 
at the origin in momentum space, whereas these signs 
are opposite for the first radial excitation. This pattern 
repeats for higher even- and odd-numbered radial excita- 
tions. Here, a direct solution of the inhomogeneous BSE 
is instructive because this equation does determine signs. 
For example, consider the pseudoscalar vertex: Fig. 6 of 
Ref. [24| illustrates a case in which the residue associated 
with the pscudoscalar meson ground-state is positive and 
that connected with the first radial excitation is negative, 
which is the behaviour found herein for lu > 0.5 GeV. The 
residue is a product of the pseudoscalar-meson's bound- 
state Bethe-Salpeter amplitude at k 2 = 0, r -(0;P 2 ), 
and p - . The latter is the expression in quantum field 
theory for the value of the Bethe-Salpeter wave func- 
tion at the origin in configuration space. Thus, the pat- 



tern exposed by the inhomogeneous BSE parallels that 
in quantum mechanics. 

It is straightforward to see that this pattern is realised 
in the second, third and fourth panels of Fig. 01 which 
depict results obtained with u> > 0.5 GeV. Therein, the 
k 2 = values of the leading amplitudes' lowest Cheby- 
shev projections are positive; and whilst that for the 
ground-state remains positive, that for the first radial 
excitation changes sign, so that it is a negative-definite 
function for k 2 > 1 GeV 2 . In performing a Fourier trans- 
form, largc-fc 2 maps onto small x 2 and hence this be- 
haviour guarantees that the Bethe-Salpeter wave func- 
tion for the first radial excitation is negative at the origin 
in configuration space. 

These observations reemphasise the peculiar character 
of the u = 0.4 GeV solution in the top-left panel of the 
Fig.[4]and explain our choice of sign for all Bethe-Salpeter 
amplitudes. The ground-state amplitude is positive at 
large-fc 2 , the first radial excitation is negative at largc- 
fc 2 , and so on. With this convention, one necessarily finds 
/4 o > 0, p^ < 0, etc., and hence, from Eq. (|2"2"]). /„.„ > 0, 
/ttj < 0. We depict the w-dependence of the leptonic 
decay constants in Fig. [5] 

The bottom four panels of Fig. 2] display low-order mo- 
ments of the exotic-pseudoscalar-meson's leading invari- 
ant amplitude, contrasted with the ground-state's zeroth 
moment. So long as ui > 0.5 GeV, the first moment of 
the exotic amplitude is bounded above by "E^o and the 
third moment is negative definite. This is the first time 
these features have been exposed but we expect them 
to be characteristic of the rainbow-ladder truncation. It 
will be important to learn whether this pattern persists 
beyond rainbow-ladder truncation. 

The top four panels in Fig. [5] compare the ampli- 
tudes of the ground-state and first-radially-excited vec- 
tor mesons. The ground-state is insensitive to lu so long 
as lu > 0.5 GeV but again the radial excitation reacts 
strongly to variations in lu. In this case, natural be- 
haviour for the excited state's amplitudes is only ob- 
tained for lu > 0.6 GeV. For smaller values, the zeroth 
moment is negative-definite and the second moment ex- 
hibits a zero. NB. The sign of the amplitudes is fixed via 
the same prescription used for pseudoscalar mesons, and 
hence f Po > 0, f Pl < 0. 

The bottom four panels of Fig. [5] display low-order mo- 
ments of the exotic-vector-meson's leading invariant am- 
plitude, contrasted with the ground-state's zeroth mo- 
ment. In this case, so long as lu > 0.6 GeV, the first 
moment of the exotic amplitude is bounded above by 
°E po and the third moment is negative definite. The sim- 
ilarity to the lower panels of Fig. 2] encourages us in the 
expectation that these features are characteristic of the 
rainbow-ladder truncation. Moreover, they suggest again 
that there is too much similarity between natural and ex- 
otic C-parity states in rainbow-ladder truncation. 

In Fig. [5] we depict the w-dependence of pscudoscalar- 
and vector-meson leptonic decay constants. Those for the 
ground-states are positive whilst those for the first radial 
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FIG. 5. Vector mesons. cj-dependence of low-order 
Chebyshev-projections of leading invariant amplitude for 
ground-, radially-excited- and exotic-states: upper four pan- 
els, ground and radial; lower four panels, ground and exotic. 
In all panels, solid - zeroth moment, ground-state; dashed 
- leading moment, comparison state; dash-dot - subleading 
moment, comparison state. Row-1, left, uj — 0.4 GeV; Row-1, 
right, uj = 0.5 GeV; Row-2, left, uj = 0.6 GeV; and Row 2, 
right, uj = 0.7 GeV. This pattern is repeated in the next two 
rows. The normalisation is chosen such that °E P0 (p 2 = 0) = 1; 
and Duj = (1.1 GeV) 3 . 



excitations are negative. The origin of this outcome in 
an internally consistent treatment of bound-states was 
explained above. Notable, too, is the small magnitude of 
the decay constant for the pion's first radial excitation: 
/ 7ri w — 1 MeV. This was predicted in Ref. [2(| and is a 
consequence of the axial- vector Ward-Takahashi identity. 
It is consistent with data on r — s- 7r(1300)^ T [6l| and 
numerical simulations of lattice- regularised QCD [231 ]. 




0.40 



0.45 



0.50 
w [GeV] 



0.55 



0.60 



FIG. 6. w-dependence of leptonic decay constants for pseu- 
doscalar and vector mesons: ground-state pion, solid line; 
radially-excited pion, dashed line; ground-state rho-meson, 
dotted line; and radially-excited rho-meson, dash-dot line. 
(Du = (1.1 GeV) 3 .) 



IV. CONCLUSION 

Using an interaction kernel that is consonant with 
modern DSE- and lattice-QCD results, we employed a 
rainbow-ladder truncation of QCD's Dyson-Schwinger 
equations in an analysis of ground-state, radially-excited 
and exotic scalar-, vector- and flavourcd-pscudoscalar- 
mesons. We confirmed that rainbow-ladder truncation is 
incapable of providing realistic predictions for the masses 
of excited- and exotic-states; e.g., the ordering between 
pseudoscalar and vector radially-excited states is incor- 
rect, and computed masses for exotic states are too low 
in comparison with other estimates. Indeed, in rainbow- 
ladder truncation, it appears that exotic states are in 
most respects too much like their C-parity partners. 

On the other hand, rainbow-ladder results do pro- 
vide information that is useful in proceeding beyond 
this leading-order. For example, in each channel the 
rainbow-ladder truncation indicates those invariant am- 
plitudes which arc likely to dominate in any solution of 
the Bethe-Salpeter equation. This knowledge can be used 
in developing integral projection techniques that sup- 
press ground-state contamination when searching for ex- 
cited states. Moreover, the response of observables, and 
the Bcthc-Salpctcr amplitudes which produce them, to 
changes in the infrared evolution of the interaction ker- 
nel can be used effectively to demarcate the domain of 
physically allowed possibilities for that evolution. This is 
valuable in qualitatively constraining the long-range be- 
haviour of QCD's /3-function. In addition, the symmetry- 
preserving character of the rainbow-ladder truncation 
and the ready access it provides to Bcthe-Salpetcr ampli- 
tudes for bound-states enable one to highlight and illus- 
trate features of hadron observables that do not depend 
on details of the dynamics. 

There arc many indications that dynamical chiral 
symmetry breaking (DCSB), of which the momentum- 
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dependence of the drcssed-quark mass-function is a strik- 
ing signal, has an enormous impact on hadron proper- 
ties. This study is one of a growing body which indicates 
that the veracious expression of DCSB in the bound-state 
problem is essential if one is to reliably predict and under- 
stand the spectrum and properties of excited and exotic 
hadrons. Achieving this will provide the power to use ex- 
tant and forthcoming data as a tool with which to chart 



the nonperturbative evolution of QCD's /3-function. 
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